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Abstract 

Cell  behavior  is  numerically  investigated  for  fast  charging  of  a  lead-acid  cell  under  constant  current  or  pulsed  current.  A  model  used 
in  the  simulation  has  been  developed  to  predict  the  cell  performance  by  considering  the  solid-state  reaction  at  the  negative  electrode.  By 
incorporating  a  rest  period  and  a  depolarization  pulse  on  high  rate  charge,  it  is  possible  to  avoid  deterioration  of  the  cell  performance  due 
to  an  increase  in  internal  resistance.  Thus,  fast  charging  of  a  lead-acid  cell  can  be  achieved  without  a  loss  of  cycle-life,  despite  the  fact 
that  higher  currents  are  forced  into  the  cell.  ©2000  Elsevier  Science  S.A.  All  rights  reserved. 
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1.  Introduction 

The  fast  charging  of  a  lead- acid  battery,  or  indeed 
other  secondary  rechargeable  batteries,  is  a  key  technology 
for  electric  vehicles.  Considerable  researches  have  been 
performed  to  reduce  the  time  required  to  fully  charge  a 
battery  without  loss  of  capacity  or  life.  To  charge  at  high 
efficiency  and  for  good  cycle-life,  it  is  important  that  the 
internal  resistance  of  the  cell  must  be  kept  low. 

Gu  et  al.  [1]  examined  the  voltage  behaviour  and  the 
variation  of  electrolyte  concentration  for  constant-current 
charge  at  25°C  and  —  18°C.  Maja  et  al.  [2]  investigated 
simulation  of  the  discharge  and  charge  processes  by  study¬ 
ing  the  effect  of  various  parameters  such  as  plate  dimen¬ 
sions,  amount  of  acid,  and  porosity  of  the  active  material. 
They  also  observed  cell  voltage  as  a  function  of  the 
charging  current  intensity  and  the  charging  modality,  e.g., 
single  or  stepped.  The  battery  was  charged  up  to  70%  of 
the  nominal  capacity,  and  then  recharged  after  a  rest 
period. 

Calasanzio  et  al.  [3]  simulated  not  only  constant-current 
charge,  for  which  the  porosity  of  the  active  material  was 
found  to  influence  the  charge  returned,  but  also  step 
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charging  with  and  without  rest  periods.  Chang  et  al.  [4] 
examined  the  temperature  rise  in  a  flooded  hybrid  battery 
during  charge.  These  authors  found  that  in  the  first  part  of 
a  rapid  charge  to  about  40%  return,  resistive  heating  was 
the  major  contributor  to  the  initial  rise  in  temperature. 
With  further  charging,  heat  from  polarization  becomes 
more  important.  Gu  et  al.  [5]  investigated  the  discharge 
and  charge  processes  by  considering  not  only  electrochem¬ 
ical  kinetics  and  mass  transport,  but  also  free  convection 
due  to  acid  stratification. 

Lam  et  al.  [6]  compared  a  pulsed-current  technique  with 
a  constant-current  recharge  for  the  rapid  charging  of 
lead-acid  cells.  They  examined  the  electrochemical  char¬ 
acteristics  and  cycle-life  performance  with  each  charging 
method.  Podrazhansky  and  Popp  [7]  used  a  depolarization 
pulse  to  reduce  the  increase  in  the  internal  resistance  due 
to  severe  concentration  polarization.  The  depolarization 
pulse  also  moderated  the  temperature  rise  of  the  cell. 
Catherino  et  al.  [8]  represented  the  charging  curves  at 
constant  current  over  a  range  of  temperature  by  combining 
the  individual  effects  of  battery  charging  and  gassing.  Kim 
and  Hong  [9]  observed  the  effect  of  the  electrode  parame¬ 
ters,  such  as  the  charge-reaction  concentration  exponent, 
morphology  parameters  and  limiting  current  density  on 
charge  performance. 

In  this  work,  a  simulation  is  made  of  the  effects  of  rest 
period  and  a  depolarization  pulse  on  the  charge  behaviour 
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Table  1 

Governing  equations  for  each  region 


Porosity  variation 
Ohm’s  law  in  solid 
Ohm’s  law  in  liquid 
Mass  balance  of  electrolyte 


Electrode  kinetics  (charge) 


(Depolarization) 
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Table  2 

Boundary  conditions  for  each  region 
Centers  of  positive  and  negative  electrode 


Vc  =  0 
Vc 

V<£2  =  0 
*2  =  0 

4>i  =  0  (positive  electrode) 
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(Negative  electrode  for  charge) 


Interface  between  positive 
electrode  and  reservoir 

e  *  Vc|+=  Vc|res 
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(Negative  electrode  for  depolarization) 


Interface  between  reservoir  and  separator 


D  •  Vc  -  CV*  |res  =  D*  ■  Vc  -  CU  *  |sep 
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Interface  between  separator 
and  negative  electrode 
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of  a  lead-acid  cell,  using  a  mathematical  model,  which 
have  been  developed  during  the  studies  of  cell  perfor¬ 
mance,  both  during  constant-current  charge  [9]  and  during 
discharge  [10,11].  The  results  of  the  simulation  will  help  to 
determine  the  best  charging  method,  which  will  reduce  the 
time  to  full  charge  and  will  not  cause  damage  to  the 
battery. 

2.  Theoretical  background 

The  governing  equations  used  to  describe  the  porous 
electrode  are  material  balance  of  a  transport  species,  ohmic 
loss  in  liquid  and  solid  phases,  and  conservation  of  charge. 
In  a  macroscopic  treatment,  the  actual  pore  geometry  is 
disregarded.  If  the  potentials  in  the  solid  and  the  liquid 
phase  are  defined,  these  quantities  are  assumed  to  be 
continuous  functions  of  time  and  space  coordinates.  In 


effect,  the  electrode  is  treated  as  the  superposition  of  two 
continua:  one  represents  the  solution,  the  other  represents 
the  matrix  [12]. 

The  electrode  reactions  during  the  charging  of  a  lead- 
acid  cell  are  as  follows. 

Positive  electrode: 

PbS04(s)  +  2H20  ->  PbO,(s)  +  HS04“+  3H++  2e“. 

(1) 

Negative  electrode: 

PbS04(s)  +  2e“+  H+-*  Pb(s)  +  HS04“.  (2) 

Lead  sulfate  is  converted  into  lead  dioxide  (positive 
electrode)  and  lead  (negative  electrode)  during  charge. 

Note  that  the  electrochemical  reaction  can  only  take 
place  where  electrons  can  be  supplied  or  removed,  which 
means  that  this  conversion  is  not  possible  on  the  surface  of 
the  lead  sulfate,  as  lead  sulfate  does  not  conduct  electric 
current.  For  this  reason,  the  Pb2+  ions  should  be  dissolved 


Table  3 


Coefficients  and  effective  properties  used  in  model  equations 
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and  transported  by  diffusion  to  the  conductive  electrode 
surface  [13]. 

The  kinetics  of  a  porous  electrode  are  affected  by  the 
local  activity,  overpotential,  electrolyte  concentration,  and 
current  density  as  a  function  of  electrode  structure.  For  a 
dissolution-precipitation  mechanism,  it  is  assumed  that  the 
electrode  reactions  in  the  negative  electrode  consist  of 
three  steps,  i.e.,  dissolution  of  lead  sulfate,  diffusion  of 
lead  ions  to  the  active  sites,  and  precipitation  of  lead  [14]. 
Thus,  the  electrode  reaction  rate  at  the  negative  electrode, 
including  the  solid-state  reaction  rate  and  the  effect  of 
electrolyte  concentration,  is  expressed  as  a  transfer-current 
density  [12].  The  transfer-current  density  at  the  positive 
electrode  can  be  obtained  from  the  Butler- Volmer  equa¬ 
tion. 


Since  a  material  balance  on  the  solid  phase  shows  how 
the  porosity  changes  with  the  extent  of  reaction,  the  varia¬ 
tion  of  the  electrode  porosity  at  each  location  can  be 
equated  to  the  density  differences  between  the  solid  prod¬ 
ucts  and  reactants  that  result  from  the  local  transfer  current 
[12]. 

With  the  assumption  of  uniformity  and  homogeneity  of 
the  electrode  in  other  directions,  the  differential  material 
balance  of  a  dissolved  species  for  a  one-dimensional  sys¬ 
tem  can  be  expressed.  Diffusivity  and  conductivity  used 
here  are  related  to  both  the  porosity  and  tortuosity  of  the 
electrode. 

The  mass  transfer  in  the  electrolyte  solution  is  given  by 
convection,  diffusion  and  migration  of  mobile  ionic  species. 
The  transference  number  and  the  partial  molar  volume  of 


Table  4 

Parameters  used  in  the  calculations 
Parameter 

Positive  electrode 
Half  thickness  of  plate 
Maximum  charge  state 
Reaction  parameter  (<2maxii0iiref) 
“at 
ac  i 

yl 

n 

Lead  dioxide  conductivity 

exl 

exml 

Negative  electrode 
Half  thickness  of  plate 
Maximum  charge  state 
Reaction  parameter  ( einVdx4i  oiref) 

aa4 
ac  4 

y4 

£4 

Lead  conductivity 

ex4 

exm4 

7lim 


Reservoir 

Thickness  of  reservoir 
Separator 

Thickness  of  separator 

Porosity 

ex3 

Electrolyte 
Acid  concentration 
Transference  number 
Partial  molar  volume  of  acid 

Diffusion  coefficient 


Conductivity 


Value 


0.0875  cm 
2620  C  cm-3 
0.073  A  cm-3 
1.15 
0.85 


0.01  (charge) 


1.0 

500  S  cm"1 

1.5 

0.5 


0.07  cm 
3120  C  cm"3 
0.1 1  A  cm"3 


1.55 

0.45 

0.01 

1.0 

4.8  X  104  S  cm"1 

1.5 

0.5 

—  102  A  cm" 3  (charge) 

—  105  A  cm" 3  (depolarization) 


0.07  cm 


0.022  cm 

0.60 

1.50 


4.9  X  10“3  mol  cm"3 
0.72 

45  mol  cm"3 

/  2174.0  2174.0 

I)  =  exp - 

\  298.15  T 


X  (1.75  +  260.0c)  X  10" 


k  =  c  X  exp 


3916.95  -  99406.0c 


721860 


1.1104  +  (199.475  -  16097.781c)c  + 
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(a) 


Fig.  1.  Charging  modalities:  (a)  constant  current;  (b)  pulse  current 
without  depolarization  current;  (c)  pulse  current  with  depolarization 
current. 


the  electrolyte  are  assumed  to  be  constant.  The  movement 
of  electrons  on  the  solid  phase  is  followed  by  Ohm’s  law. 
The  thickness  of  the  electric  double  layer  is  smaller  than 
the  pore  volume,  so  it  is  assumed  to  be  electrically  neutral 
for  the  electrolyte  solution  [15].  The  divergence  of  total 
current  density,  therefore,  is  zero.  All  the  charge  leaving 
the  solid  phase  enters  the  solution  phase. 

The  polarization  equation  is  used  to  describe  the  voltage 
jump  at  the  interface  between  the  solid  and  liquid  phase 
and  the  dependence  of  the  electrolyte  concentration  on  the 
local  reaction  rate.  In  the  reservoir  and  the  separator,  the 
transfer-current  density  is  set  at  zero.  The  governing  equa¬ 
tions  used  in  the  numerical  simulation  are  summarized  in 
Table  1  [9-11]. 


The  initial  values  for  electrolyte  concentration  and 
porosity  are  obtained  from  the  results  of  a  discharged  cell 
before  the  start  of  the  charging  process.  The  initial  poten¬ 
tial  distribution  is  calculated  from  the  kinetic  equation  of 
electrode  reaction. 

At  the  current-collectors,  as  the  battery  cell  is  a  closed 
system,  the  potential  and  the  electrolyte  concentration  are 
set  to  constant  values  [5].  The  set  of  boundary  conditions  is 
given  for  the  molar  flux  of  the  mobile  ionic  species,  charge 
conservation,  and  porosity  change.  The  boundary  condi¬ 
tions  used  to  solve  the  governing  equations  are  summa¬ 
rized  in  Table  2,  and  the  coefficients  used  in  the  develop¬ 
ment  of  the  model  are  presented  in  Table  3. 

The  equations  given  in  Tables  1  and  2  are  solved  by 
numerical  methods.  The  space  derivatives  are  discretized 
into  the  method  of  finite  differences  by  means  of  three-point 
approximations,  and  the  time  derivatives  are  formulated  by 
the  Crank-Nicolson  method.  The  set  of  non-linear  equa¬ 
tions  in  the  multi-region  system  is  solved  by  the  iterative 
Newton  method  and  modified  Newman’s  MB  AND  [16]. 
The  material  properties  and  the  cell  parameters  used  in 
simulation  are  presented  in  Table  4  [9], 


3.  Results  and  discussion 

The  model  used  in  the  simulation  includes  the  dissolu¬ 
tion-precipitation  mechanism  in  the  reaction  kinetics  of 
the  negative  electrode.  The  cell  used  in  the  simulation  was 
numerically  discharged  with  1.7  mA  cm-2  for  8  h  at  25°C, 
and  then  experienced  rest  period  for  1  h  at  the  same 
temperature. 

The  three  charging  methods  used  in  the  study  are 
illustrated  in  Fig.  1.  Fig.  la  shows  that  the  cell  is  charged 
with  a  constant  current  to  the  voltage  cut-off.  The  schedule 
in  Fig.  lb  is  for  pulsed-current  charging;  current  Ic  is 
forced  into  the  cell  for  time  tc  and  then  a  rest  period  is 


Fig.  2.  Predicted  voltage  behaviour  for  constant-current  charge  at  30  mA 

-2 

cm  . 
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Fig.  3.  Predicted  profiles  of  acid  concentration  for  constant  current  at  30  Fig.  5.  Predicted  profiles  of  acid  concentration.  Charging  current  is  30 
mA  cm  2 .  mA  cm-2;  duration  of  charge  pulse  is  27  s;  rest  period  is  3  s. 


followed  for  time  tR.  A  series  of  charging  pulses  and  rest 
periods  are  repeated  until  the  cell  voltage  reaches  the 
voltage  cut-off.  The  rest  period  has  a  role  in  relaxing  the 
concentration  gradient  or  polarization  in  the  cell,  which 
reduces  the  temperature  rise  of  the  cell.  Fig.  lc  illustrates 
that  a  depolarization  current  7d  is  used  for  time  td  at  the 
end  of  the  rest  period.  Pulsed-current  charging  with  no 
depolarization  pulse  requires  a  longer  time  to  eliminate  the 
concentration  gradient.  Because  a  rest  period  of  short 
duration  has  little  effect  on  the  reduction  of  the  concentra¬ 
tion  polarization,  a  depolarization  pulse  is  used  to  decrease 
the  concentration  gradient. 

The  variation  of  cell  voltage  with  time  for  charging  at  a 
constant  current  of  30  mA  cm~2  to  2.6  V  is  shown  in  Fig. 
2.  The  slope  of  the  voltage  increment  increases  with 
progressing  charge,  and  at  the  end  of  charging,  its  value 
rises  rapidly.  This  is  because  the  internal  resistance  of  the 
cell  becomes  large.  Charge-acceptance  or  charge-current 
efficiency  describes  the  share  of  the  current,  which  is 
actually  accepted  by  the  battery  and  can  be  retrieved  by 


Fig.  4.  Effect  on  rest  period  on  cell  voltage.  Charging  current  is  30  mA 
cm-2. 


subsequent  discharging  [13].  This  value  is  not  constant,  but 
depends  largely  on  the  state  of  charge,  i.e.,  charge-accep¬ 
tance  decreases  with  progress  of  the  charge  reaction.  To 
ensure  that  the  charge-acceptance  has  a  large  value,  a 
higher  surface  area  for  the  electrode  and  a  lower  internal 
resistance  of  the  cell  is  required. 

The  concentration  distribution  for  constant-current 
charge  under  the  same  conditions  as  in  Fig.  2  is  displaced 
in  Fig.  3.  At  the  start  of  charge,  the  electrolyte  is  uniformly 
distributed  in  the  cell,  but  the  electrolyte  concentration  in 
the  positive  and  negative  electrode  increases  much  faster 
than  that  in  other  regions.  Thus,  for  a  constant-current 
charge,  the  concentration  distribution  is  not  uniform,  espe¬ 
cially  at  the  interface  between  the  positive  electrode  and 
the  reservoir,  as  well  as  at  between  the  negative  electrode 
and  the  separator.  In  these  two  regions,  the  concentration 
gradients  become  steeper. 

In  order  to  charge  a  battery  quickly,  a  higher  current 
must  be  forced  into  it.  One  of  the  methods  that  increases 
the  charge-acceptance  and  reduces  the  internal  resistance  is 


Fig.  6.  Predicted  profiles  of  acid  concentration.  Charging  current  is  30 
mA  cm-2 ;  duration  of  charge  pulse  is  27  s;  rest  period  is  6  s. 
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Fig.  7.  Predicted  profiles  of  acid  concentration.  Charging  current  is  30 
mA  cm-2 ;  duration  of  charge  pulse  is  27  s;  rest  period  is  23  s. 


pulsed-current  charge.  A  high  current  is  passed  to  the 
battery,  which  is  followed  by  a  rest  period  to  disperse  ions 
throughout  the  electrolyte.  The  effect  of  the  rest  period  on 
the  cell  voltage  with  time  is  given  in  Fig.  4  for  pulses  of 
constant  magnitude.  The  time  is  fixed,  at  27  s  and  the 
relaxation  time  is  3  and  6  s,  and  for  the  changing  pulse  23 
s.  In  effect,  the  slope  of  the  voltage  increment  is  lower 
with  a  rest  period  and  the  value  decreases  with  increasing 
duration  of  the  rest  period. 

The  concentration  distributions  for  pulsed-current  charge 
with  the  same  conditions  as  those  used  in  Fig.  4  are  given 
in  Figs.  5-7.  For  a  rest  period  of  3  s  (Fig.  5),  the 
concentration  polarization  throughout  the  electrolyte  is  re¬ 
duced  by  more  than  that  without  a  rest  period  (cf.  Fig.  2). 
When  the  duration  of  rest  period  is  6  s  (Fig.  6),  the 
concentration  gradient  between  the  negative  electrode  and 
the  separator  is  lower  than  that  in  Fig.  5,  but  that  between 
the  positive  electrode  and  the  reservoir  is  still  large.  For  a 
rest  period  of  23  s  (Fig.  7),  the  behaviour  of  the  concentra¬ 


Fig.  8.  Effect  of  depolarization  pulse  on  charge  behaviour.  Charge  current 
is  30  mA  cm”2  and  duration  is  27  s.  Rest  period  is  3  s  and  duration  of 
depolarization  pulse  is  2  s. 


Fig.  9.  Predicted  profiles  of  acid  concentration.  Charging  current  is  30 
mA  cm”2  and  duration  of  charge  pulse  is  27  s.  Rest  period  is  3  s  and 
duration  of  depolarization  pulse  is  2  s;  depolarization  current  is  —  30  mA 
cm-2. 


tion  distribution  is  similar  to  that  shown  in  Fig.  6  for  a  6-s 
rest  period. 

The  effect  of  a  depolarization  pulse  on  the  charging 
voltage  is  presented  in  Fig.  8.  During  charging,  positive 
ions  move  to  the  negative  electrode  and  negative  ions 
move  to  the  positive  electrode.  The  ions  arriving  at  the 
electrode  progressively  shield  the  electrode  from  further 
charge-transfer.  For  a  pulsed-current  charge  without  a 
depolarization  pulse,  the  rest  period  reduces  the  concentra¬ 
tion  gradient.  As  the  duration  of  the  rest  period  is  in¬ 
creased,  however,  it  takes  more  time  to  achieve  full  charge. 
Thus,  a  depolarization  current  serves  to  eliminate  the 
concentration  polarization  and  to  reduce  the  charging  time. 
Note  that  the  duration  of  depolarization  pulse  is  smaller 
than  that  of  charging  pulse.  Therefore,  the  depolarization 
pulse  should  not  cause  any  significant  discharge  of  the 
battery.  The  duration  between  the  charging  and  depolariza- 


Fig.  10.  Predicted  profiles  of  acid  concentration.  Charging  current  is  30 
mA  cm”2  and  duration  of  charge  pulse  is  27  s.  Rest  period  is  3  s  and 
duration  of  depolarization  pulse  is  2  s;  depolarization  current  is  —  60  mA 
cm-2. 
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Fig.  11.  Predicted  profiles  of  acid  concentration.  Charging  current  is  30 
mA  cm-2  and  duration  of  charge  pulse  is  27  s.  Rest  period  is  3  s  and 
duration  of  depolarization  pulse  is  2  s;  depolarization  current  is  —90  mA 
cm-2. 


tion  pulse  ( tR ,  Fig.  lc)  is  set  to  3  s  and  the  charging  pulse 
is  applied  at  30  mA  cm-2  for  27  s.  Whereas  the  depolar¬ 
ization  pulse  is  forced  to  the  cell  for  2  s,  the  current 
density  is  changed  to  —30,  —60  and  —90  mA  cm-2.  As 
the  current  density  of  the  depolarization  pulse  is  increased, 
the  voltage  increment  is  decreased.  This  is  because  the 
internal  resistance  of  the  cell  is  reduced. 

The  concentration  distributions  for  pulsed-current  charge 
with  a  depolarization  pulse  under  the  same  conditions  as  in 
Fig.  8  are  presented  in  Figs.  9-11.  In  Fig.  9,  the  depolar¬ 
ization  pulse  is  set  at  —30  mA  cm'2  and  gives  similar 
results  to  those  obtained  with  a  rest  period  of  long  dura¬ 
tion,  as  shown  in  Figs.  4  and  5.  When  the  intensity  of  the 
depolarization  pulse  is  —60  mA  cm-2  (Fig.  10),  the 
concentration  gradient  is  reduced  more  than  that  in  Fig.  9. 
The  concentration  polarization  is  relaxed  even  further  at 
—  90  mA  cm'2,  as  shown  in  Fig.  11.  Using  the  depolariza¬ 
tion  pulse,  it  can  be  seen  that  the  internal  resistance  is 
reduced  and  the  charge-acceptance  is  increased. 


4.  Conclusions 

Variations  in  cell  voltage  and  concentration  distribution 
in  a  lead-acid  cell  under  constant  or  pulsed-current  charge 
have  been  investigated  using  a  mathematical  model,  which 
includes  the  dissolution-precipitation  mechanism  in  the 
negative  electrode.  The  results  of  applying  the  model  can 
be  summarized  as  follows. 

(i)  For  a  constant-current  charge,  the  internal  resistance 
of  the  cell  has  a  greater  effect  on  the  charge  performance. 
The  cell  voltage  rises  rapidly  at  the  end  of  charging. 

(ii)  The  slope  of  the  voltage  increment  is  lower  when  a 
rest  period  is  included.  Its  value  decreases  with  increasing 
duration  of  the  rest  period. 


(iii)  Using  a  pulsed  current  with  a  depolarization  pulse, 
the  internal  resistance  is  reduced  and  the  charge-accep¬ 
tance  is  increased.  Therefore,  a  high-rate  charge  can  be 
performed  without  an  increase  in  internal  resistance  or 
reduction  in  charge-acceptance. 

5.  List  of  symbols 

a  active  surface  area  of  electrode  (cm2  cm-3) 

flj,  a2  coefficients  used  in  governing  equations  and 
boundary  conditions 

amax  maximum  specific  active  surface  area  of  elec¬ 
trode  (cm2  cm'3) 

c  concentration  of  binary  electrolyte  (mol  cm'3) 

cref  reference  concentration  of  binary  electrolyte  (mol 
cm'3) 

D  diffusion  coefficient  of  binary  electrolyte  (cm2 

s'1) 

F  Faraday’s  constant  (96,487  C  mol'1) 

i  total  applied  current  density  based  on  projected 

electrode  area  (A  cm'2) 

?!  current  density  in  solid  phase  (A  cm'2) 

i2  current  density  in  conducting  liquid  phase  (A 

cm'2) 

zolrcf  exchange  current  density  at  cref  for  positive 

electrode  (A  cm'2) 

zo4ref  exchange  current  density  at  cref  for  negative 

electrode  (A  cm'2) 

j  reaction  current  per  unit  volume  of  electrode  (A 

cm'3) 

VjIIT1  limiting  current  density  for  negative  electrode  (A 
cm'3) 

MW;  molecular  weight  of  species  i  (g  mol'1) 

Q  charge  density  in  electrode  (C  cm'3) 

<2max  theoretical  maximum  capacity  (C  cm'3) 

R  universal  gas  constant  (8.3143  J  mol'1  K'1) 

t  time  (s) 

t'l  transference  number  of  H+  with  respect  to  sol¬ 

vent  velocity 

T  absolute  temperature  (K) 

UPb0  equilibrium  potential  at  cref  for  positive  elec¬ 

trode  (V) 

v  *  volume-average  velocity  (cm  s'1) 

x  distance  from  the  center  of  positive  electrode 

(cm) 

Greek  letters 

aal  anodic  transfer  coefficient  for  positive  electrode 

otl  cathodic  transfer  coefficient  for  positive  elec¬ 

trode 

oa4  anodic  transfer  coefficient  for  negative  electrode 

acA  cathodic  transfer  coefficient  for  negative  elec¬ 

trode 

y  1  concentration  exponent  for  positive  electrode 

y4  concentration  exponent  for  negative  electrode 

s  porosity 


S.C.  Kirn,  W.H.  Hong  /  Journal  of  Power  Sources  89  (2000)  93-101 


101 


£sep  porosity  of  separator 

£pb0,,ini  porosity  of  positive  electrode  at  initial  state  of 
charge  reaction 

s pb;ini  porosity  of  negative  electrode  at  initial  state  of 
charge  reaction 

£  1  morphology  parameter  for  positive  electrode 

£4  morphology  parameter  for  negative  electrode 

17  total  local  overpotential  with  respect  to  equilib¬ 

rium  potential 

k  electrolyte  conductivity  (S  cm-1) 

p,  density  of  species  i  (g  cm-3) 

(X,  conductivity  of  electrode  matrix  (S  cm-1) 

(f)  j  potential  in  electrode  matrix  (V) 

cf) 2  potential  in  solution  (V) Superscripts 

ex  exponent  on  porosity 

exm  empirically  determined  constant  for  tortuosity  of 

the  solid  matrix 
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